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Abstract. We report 2D time-dependent non- linear magneto- hydrodynamical simulations of 
waves in the atmospheres of roAp stars. We explore a grid of simulations in a wide parameter 
, space. The aim of our study is to understand the influence of the atmosphere and the magnetic 
field on the propagation and reflection properties of magneto-acoustic waves, formation of shocks 
and node layers. 
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RoAp stars are late-A, chemically peculiar stars with effective temperatures between 
6500-8100 K and global dipolar-like magnetic fields with strengths between 1-25 kG. 
They show low-order angular degree p-mode pulsations with periods between 4 and 22 
i-G ■ minutes. These non-radial pulsations are aligned with dipolar field axis. Magneto-acoustic 
oscillations in peculiar A stars are of particular interest due to unique opportunities to 



We solve the basic 2D non -linear adiabatic equations of t h e ideal MHD by means o f the 
numerical code described bv lKhomenko fc Colladosl (l2006h : Kho menko et al. I (|2008h . We 



study the interaction of pulsations, chemical inhomogeneities, and strong magnetic fields. 
+3 [ Recen t reviews on the properties of these stars a nd their pulsations can be found in Kurtz 
(|2008h : ICunhal (E007h : iKochukhwl (|2007l I2008TI . 

assume that: (1) the magnetic field varies on spatial scales much larger than the typical 
■ wavelength, allowing the problem to be solved locally for a plane-parallel atmosphere 
with a homogeneous inclined magnetic field; (2) waves in the atmosphere are excited 
by low-degree pulsation modes with radial velocities exceeding horizontal velocities. The 
unperturbed atmospheric model has an effective temperature of T c g — 7750 K and gravi- 
tational acceleration at the surface logio.g = 4.0. The simulation grid covers the magnetic 
field strength B = 1 — 7 kG; magnetic field inclination to the local vertical 7 = — 60 
degrees and driving periods T = 6 — 13 minutes. An example of the wave pattern devel- 
oped in the simulations is given in Fig. 1. We use the longitudinal and transversal field 
projections of the velocity in order to separate clearly the fast and slow MHD waves. 
Our first results pick up some observed properties of the roAp stars pulsations, such as: 
rapid growth of the wave amplitude with height, presence of the node surfaces, large 
variations in pulsation properties depending on the parameters of the model. Velocity 
signal observed in the upper atmospheric layers of roAp stars is mostly due to running 
slow mode (acoustic) waves propagating along the inclined magnetic field lines. The node 
structures and the rapid phase variations at the lower atmospheric layers are due to mul- 
tiple reflections and interference of the slow and fast mode waves. The disc-integrated 
velocity signal produced by the atmospheric pulsations of such a star will depend in a 
complex way on the inclination of the magnetic axis with respect to the observational 
line of sight and will be a subject of our further study. 
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Figure 1. Height-time variations of the longitudinal (parallel to the field, left) and transversal 
(perpendicular to the field, right) velocities for B=l kG, T=360 sec at latitudes where the 
inclination equals 0, 30, 45 and 60 degrees. The color bars give velocity scale in km/sec. Zero 
height corresponds to the photospheric base. Dotted lines marked with numbers are contours 
of constant c| /v\ . Height dependences of the sound speed cs (solid line) and Alfven speed va 
(dashed line) are plotted over the top panel, the scale is given by the upper axis. 
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